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\o 

We propose the description of the granular matter which is based on distribution of dry friction 
coefficients. Using such a concept and a simple one-dimensional packing of grains we solve the silo 
C^l , problem. The friction coefficients at contacts are determined both by geometry of packing config- 

» ! ■ uration and the stress distribution in a medium. Within such an approach the Janssen coefficient 

kj is determined and its dependence on the particle-particle and boundary-particle friction coeffi- 
cients is obtained. Also we investigate the conditions for the appearance of the maximum in the 
pressure distribution with the depth with overweight on top. As an outcome of our work we propose 
the general framework to the description of the granular matter as a continual medium which is 
£Nj ' characterized by the field of the dry friction tensor. 
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Introduction 



The properties of granular matter differ drastically from those of other continuous media like solids, liquids and gases. 
I | Despite its mechanical nature the problem of description of granular media is still open problem. Some arguments 
put forward recently Q, Q even question the possibility of the description of such media basing on the hydrodynamic 
approach. The latter usually applied to continuous media with short-range interparticle forces. But nonconservative 
nature of the friction force hampers direct application of such an approach based on local balance equations for 
physical quantities like momentum, energy, entropy, etc. In such a complex situation with the dynamics the static 
properties of the granular medium are easier to investigate. There are several characteristic static phenomena, which 
any theory should explain. They are connected with the static distribution of the pressure along a silo. Namely, the 
deviation from Pascal law and the nonmonotonic dependence of the apparent mass on the overweight on the top of a 
silo. Lets us review them shortly. 
' In contrast with conventional fluids, "hydrostatic" pressure in granular media reaches the asymptotic finite value 
£T) . po at some finite depth A, and can be described as 

O ' 

p(*)=Po(l-exp(-*/A)), A = ^j (1) 



where R is the width of the silo, \i is the static friction coefficient between the grains and the walls of the silo, and 
kj is the so-called Janssen coefficient. First explanation of this behavior originated from the Jannsen model Q (see 
also jjj]), based on two simple assumptions: i) the linear relation between cr xz and a xx components of the stress 
tensor, which is analogous to the Coulomb- Amonton law for the dry friction, and ii) the linear ratio between other 
two components of the stress tensor: a xx — kj<j zz . Since then, many approaches to the statics of granular medium 



fj ' have been proposed either giving the grounds for the Janssen relations |a, l| , or building the theory without using 



them 

One of the simplest approach treats the granular media within the framework of linear isotropic elasticity. Indeed, 
the Janssen coefficient kj might be expressed through the Poisson ratio of the media, and numerics recovers the 
Janssen relations between the components of the stress tensor [fj]. Detailed comparisons with experimental results 
even allow to extract the relation between the Poisson ratio and the granular packing fraction [(|. However, such a 
picture is hardly adequate for the granular media of rigid particles. In addition, it is hard to control realization of the 
Coulomb threshold condition everywhere at the walls, which, obviously, influences the experimental results 10] . 

Other approaches do not exploit elastic nature of the grains, but use the Coulomb- Amonton law for the stress tensor 
components inside the media or at the walls |ll| . To compete with indeterminancy of the problem, the granular media 
is assumed to be at the verge of Coulombic failure everywhere in the bulk. However, this assumption about fully 
mobilized friction seems to be too restrictive 4], since the frictional forces can be varied provided that the medium 
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stays at rest. In order to simplify the problem which include randomness in the distribution of such coefficients, 
specific assumptions about the geometrical distribution of the grains should be made, e.g. within some lattice models 
d, 0] • Such models give the results which are consistent with continuum theories for the average stresses. 

A particularly interesting issue is the weight distribution in a silo with some overweight on top. Experiments show 
the maximum in the pressure distribution with the depth, in contradiction with predictions of the simple Janssen 
model 11]. Interestingly, this maximum appears almost at the same depth where pressure at the silo without the 
overweight changes in e times. Although some of mentioned approaches do not contradict experiment, but the physical 
reasons of this phenomenon are still unclear. 

The aim of this paper is to propose the description of the granular media which incorporates an additional dry 
friction tensor field. In the simplest ID geometry considered below it reduces to the distribution of the dry friction 
coefficients. We show that in this case it is possible to get the expression for the Janssen coefficient. Besides, we arc 
able to explain the the weight distributions in experiments with overweight as a result of the inhomogeneity in the 
friction coefficient distribution. The existence of characteristic scale for such inhomogeneity leads to the parametric 
dependence of the relation "apparent mass - filled mass " on the ratio of two length scales characteristic inhomogeneity 
and the saturation length for the bulk pressure. 

The structure of the paper is as follows. In Section [I] we introduce the /x-model and illustrate it on the simple ID 
granular packing, which resembles the real granular packing in silo. The results and comparisons are presented in 
Section ^of the paper. We discuss the force distribution with and without overweight for simple packing and also 
find how Janssen coefficient kj can be evaluated from microscopic characteristics of the media. In Section IIIII we 
discuss the grounds of /x- model using the results in Section [D] Conclusions are given in Section llVl 



I. /i-MODEL FOR SIMPLE ID PACKING 

The grains at contact are subjected to static friction force. The value of such a force is determined both by the 
geometry of the contact surfaces and the stress at the contact. The value of the friction force is between and the 
maximum value. The last is described by the Coulomb- Amonton law for static friction. Exact value is defined by the 
condition to keep the grain at rest. Since a grain is in contact with its neighbors at several points the net force is 
the sum of reaction forces for each contact. Because of this, granular system cannot be described using deterministic 
approach. There are more unknown variables than the equations. The granular matter can implement different force 
configurations even if the packing configuration is the same. 

Between two grains the reaction force acts. It is expedient to split it into two components: normal (normal reaction 
force) and tangential (friction forces). For each point of contact we can write 

F T = n-F n , (2) 

where coefficient /i is different for different contact point but it's value is in the interval [0,/x/], where \if is the 
static friction coefficient and determines the static friction angle. Obviously these coefficients for each contact point 
together with the normal stress determine the force configuration of the system. In the continuous limit, the set of 
these coefficients transforms into the tensor field /j,. Such a tensor field becomes additional characteristic of a granular 
medium. 

To illustrate the idea we apply this model to simple granular packing shown in Fig. ^ It is a regular pseudo-ID 
packing of identical rigid spheres in the vertical chute with a proper width. We choose 1-D model in order get explicit 
analytical results. Though such models seem to be oversimplified they grasp the main feature of the granular systems, 
namely the jamming. Also they are widely used in modelling avalanches within framework of SOC 12]. In particular 
as is shown in 1] 1-D systems shows nontrivial dynamic which is due to nonpotential character of the interparticle 
interactions. Moreover, according to the coefficient of friction [i is defined in a way irrespective on the dimension 
of the system. Thus it can be defined both for the grains and for the mesoscopic regions of the granular matter thus 
giving rise to the introduction of the /i-field. 

Denote the normal force between grains i and i + 1 by P^j+i, and between the grain i and the wall by iVj. Let us 
denote the friction coefficients for each "grain-grain" pair by /i»,t+i, and the ones for the "grain- wall'" pair by <x,-. The 
set of values -Pj,i+i, Ni, fJ,i,i+i, and cr, determines the force configuration of the system. The geometrical configuration 
of the system is represented by the angles between normal to the surface of contact of two grains and the vertical as 
Oj^+i. The force equilibrium for each grain can be written as 



Ni - I', :.,>.>, ,., - Pi,i+i$i,i+i = 

-rrug + atN - /', ..,(), i., + P, 4 +i6» M+ i = 



(3) 



3 



where 



FIG. 1: Pseudo one-dimensional packing of identical spheres 



= sin - fi^i+i COSa,,, + i , 

Oi,i+l = cosa^i+i + (H,i+i sina^i+i 



We omit the equation for the moment of force to illustrate how /i-approach works. The ID nature of the system 
permits to write the expression for the the force between two particles, namely: 



^-y: mkTk+u +iwb. 



1.0 



$i,i+l f~^ Vk,k+1 + (7k 



1?i,i+l 



(4) 



where 



n 



<in—\,n ®n 
Vrhn+1 + Q~n 



and 



1 + Hi,i+itgaj,i+i 

tg&i,i+l — A'i.i+l 



and Po is the overweight. 

Equation (0} allows to analyze the load distribution through the silo with and without overweight with some 
distribution of the contact friction coefficients 



II. RESULTS OF THE /i-MODEL FOR SIMPLE ID PACKING 



It is well known that force distribution in granular media depend not only on the silo height but also on the history 
and method of preparation of the granular matter sample |13L Il4j . Within the approach proposed such a method can 
be modelled, e.g. by the distribution of the coefficients <7j, which shows the stresses at the walls. Basing on the Eq. 
(QJ as the exact microscopic solution of the model problem, we can model such a feature, by the set of coefficients 
{/^M+i: a i} an d angles {a i}i+ i}. These data determine the geometrical and stress configuration of the system. 

The question about deviation of static pressure distribution in granular media from the Pascal law is of particular 
interest since it is the characteristic feature of such kind of materials. Note that that Eq. (0J reduces to the Pascal law 
with the pressure being proportional to the depth if the walls are absolutely smooth, i.e. <Ji = 0. This limiting case 
is in correspondence with the fact that the deviation from the PL is due to nonlinear dependence of the tangential 
component of the stress along the wall. 

Other configurations give the deviation form the PL. We investigate some of them numerically. To compare our 
results with known experimental data we take the random distributions of the friction coefficients Oi £ [0.25, 0.27] and 
Mm+i [0.63, 0.73] with the limiting values chosen as the best fits to the simulational values of 6], and ctj j+j = 7r/4. 
Statistical averaging was performed over the 100 configurations for 50 grains, each one of mass 10 g. The resulting 
dependence between apparent and filled mass is shown in Fig. [21 



4 




FIG. 2: Dependence the apparent mass on the filled mass of Eq. (|1J with random distributions of £ [0.63, 0.73] and 

Oi G [0.25, 0.27], and a^i+i = it/4, mi = Wg. Data are comparable with results of Ref. [g. 



A. Force distribution in the silo without an overweight in continuous limit 

Another relatively simple case is the "homogeneous" granular packing with fj-ij+i = fi, <Ji = \x w and a^i+i = a. In 
this case one can take the continuous limit in Eq. For a silo without an overweight (Pq — 0), one gets 

P(h) = pg\ (l - e-^) , (5) 

with characteristic length 

A= j— (6) 

where h is the effective height: h = h ■ $ ■ (fi w + rf), and configuration parameters 

1 + fi ■ tan a 1 + sin a 

7/ = , u = sin a — [i cos a , Q — 

tan a — [i cos a 

and d is chute size. 

So, we can see how the Janssen result can be obtained due to microscopic approach. 

Note, that if tana = /x, Eq. (0 again reduces to the PL since grains do not lean against the wall. Another words 
the static friction between grains is enough to keep them at rest without any wall. 

P{h) = pgh ■ cos a . (7) 

Our result can be compared with Janssen formula Formula Q was obtained for vertical cylindrical chute. To 
compare it with results of our model it should be obtained for parallelepiped chute, which is infinite in one direction 
and has the profile as shown on Fig. ^ Calculation for this case transforms formula to the following: 

p(z) = po (1 - exp(-z/A)) , A = -4r (8) 

In the "homogeneous" state the analogue of the Janssen coefficient kj can be devised via comparison of Eqs. © 
and 10 with the result JSJ: 

^,, /M „ + 1) ,„(i^±i). (9) 

Equation 10 relates the Janssen coefficient with the "microscopic" characteristics of the granular packing. In Fig. |3| 
we illustrate the dependence of kj on friction at the wall (fJ, w ) for different values of internal friction \i. 
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FIG. 3: The Janssen coefficient kj as a function of friction [i w on the walls from Eq. @ for different \i. 

B. Force distribution in the silo with overweight 

The influence of the /i distribution on the distribution of the pressure can also be shown by considering the system 
when overweight is present. 

For our numerical studies of the force distribution in the silo with overweight we use Ci = 0.25 and value of 
overweight Pq = 80. 8g, which correspond to the value of grain- wall friction coefficient and overweight of Ref. 
Vessel contains 50 grains, each one of mass 10 g. The resulting dependence between apparent and filled mass is 
shown in Fig.0] Bottom curve shows the force distribution in "homogeneous" media, where all Hi.i+i = Ho without 
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FIG. 4: Dependence of the apparent mass on the filled mass. Bottom curve: force distribution without overweight. Middle 
curve: force distribution with 80.8 g overweight for /ii, j+i = 0.65. Top curve: force distribution with 80.8 g overweight for 
Mi, s+i = 0.65 — O.le - 0,5 \ For all curves <7; = 0.25, a;,i+i = tt/4, mi = W(g). Data are comparable with results of Ref. |g. 

overweight. As one can see from previous section this curve is nothing but Janssen exponential distribution. In 
case of overweight presence in such a " homogeneous" media, the force distribution is described by the middle curve. 
This result was also predicted by Janssen but it contradicts with the experiment. Since the force configuration is 
determined by the distribution of the coefficients m^+i, we suppose that the latter has the same functional character 
as that for the pressure without overweight. The grounds of this assumption will be expanded in the next section. 
Choosing the ^-distribution as following: 

MM+i =Vo + be- c \ (10) 

we adjust parameters b, c so that to achieve the best fit (the top curve on Fig.^J to data of Ref. 6]. Thus we can see, 
that the maximum in the force distribution in the silo with overweight can be obtained if the friction coefficients are 
changed in the same way as the pressure changes without an overweight. 
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III. DISCUSSION OF THE /i-MODEL 

The results obtained in previous section for simplified ID model within the framework of //-model can be extended 
further since the final results do not contain any microscopic characteristics of such a specific model. Indeed, the 
proposed approach allows us to switch the description of the granular media from the consideration of force network 
to the distribution of friction coefficients, which in continuous limit transforms to the field of tensor fi. Such a field 
becomes an additional characteristic of granular media which cannot be obtained due to common approach (e.g from 
Newton's equations). Additional statistical arguments about how this field is distributed should be used. 

As one can see from Sec. |D] we modelled the distribution of the u in two different ways. First was uniform 
distribution /i = const, and we showed how such an assumption conforms with previous theoretical results and 
experimental data. Another distribution was of exponential form and here we give the grounds for such a choice. 

A. Spatial distribution of the /t-field 

Let us consider the distribution of pressure in a silo with additional weight on its top. Suppose, that overweight is 
implemented by another silo of the same material, in which the pressure has reached its saturated value. Thus, the 
considered part of the granular media is actually belongs to the region where the pressure has reached it's saturated 
value. Therefore the pressure in this part must be equal to overweight. These reasonings are confirmed by Janssen's 
model 0] and our results, in which it is assumed that friction either at the wall or in the bulk is const. This is in 
obvious contradiction with the experiment which shows that there is a maximum in the pressure dependence on 
depth. 

It is possible to get such a nonmonotonic dependence of pressure on overweight in some theoretical approaches (see 
e -g- 0)i but there is no clear explanation why it appears. In addition, they are based on the assumptions fi = const, 
which is adequate only near the Coulomb threshold. As we can see from the experiment, there are the "screening" 
region of size A. In case of absence of the overweight, the pressure increases there and in case of overweight presence, 
the pressure has maximum within this region. Under overloading the stress configuration within this region changes 
more drastically than in the bulk. Within the approach proposed it can be described as the spatial distribution of 
the coefficients and o~i, in general. In continuum approach it corresponds to some dependence fi{z) in the bulk 

and o~{z) at the boundary. 

Since fi characterizes the stress configuration, its spatial character must be similar to that of pressure in the system 
without overweight. This assumption can be viewed as a first approximation in the expansion spatially distributed 
function in a series of approximation. Indeed the function fi{z) can be written as 

fi(z) = no + n^z) + fi 2 (z) + . . . . (11) 

Here /to is constant part of fi, which is used in most of theoretical descriptions. Within the approach proposed fii(z) 
changes in the same way as the pressure without overweight changes. Then fi 2 takes into account changing in pressure 
in granular medium where fi(z) = /io + fii(z), etc. 

As one can see from previous section, to satisfy the experimental results, it is enough to model the distribution of 

m,i+l as: 

IM,i+i = fJ>o + be~ a , (12) 

omitting higher order approximations. We take /io as the bulk value, coefficient b itself governs the gap between bulk 
and border value of fi, both b and c are responsible for the speed of /i increasing. At the values of parameters, found in 
Sec.m fi becomes saturated very quickly, so we have very thin region where fi changes from 0.55 to 0.65 (see Fig.|S}. 

But, as one can see from Fig this region, where fi has such a dependence is enough to change the pressure 
distribution from the flat curve to the curve with the maximum. There is no maximum in the pressure dependence 
with depth, when fi = const. This result, which was predicted by Janssen model, is confirmed by physical arguments, 
given above. Maximum appears when fi is distributed within a bulk in a given way. To satisfy the experimental 
results, region where fi changes in such a way must be narrow. Thus we can say about changing only near-boundary 
values of fi. 

Note, that we have an essential difference in pressure dependencies only if overweight is present. When there is no 
overweight, the pressure dependence for constant fi and the pressure dependence for changing fi are almost the same. 
We can see that presence of an overweight reveals the real distribution of fi. 

So we can conclude the following: a) the value of fi is not constant and has some distribution; b) such a distribution 
can be revealed only by the presence of an overweight and this distribution is essential only near the boundary. 
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FIG. 5: Dependence of the spatially distributed n on the depth in units of filled mass. 

The last fact can be explained in the following way. When we fill the silo with the granular material some stress 
configuration is implemented in it. If then we put an overweight at the top of the silo, that configuration will become 
broken, and new one will be implemented. This happens because near-boundary layer feels this overweight and react 
on it. Other layers do not feel the overweight in essential way because of jamming in the upper layer. Other words, 
to fill the silo to height 2h at once is not the same that to fill the silo to height h first and then h again. This also 
illustrates how the stress distribution in the granular media depends on the history of packing creation. 

B. Macroscopic parameters 

Previous consideration allows to conclude that there must be at least two characteristic scales which characterize 
the state of silo under overloading. Followed by [|| we built rescaled dependencies (see Fig. EJ). 




12 3 4 



M ,i„'M,„ 

FIG. 6: Rescaled dependencies the apparent mass on the filled mass for different parameters b, which changes from —0.2 to 
through 0.02. The bottom curve corresponds to b = 0, the upper one corresponds to b = —0.2. Overweight and saturated 
pressure equals 80.32 g. Apparent and filled mass expressed in overweight units. Other parameters are the same as for Fig. [I] 

As one can see the bottom curve corresponds to the fj, = const because the fact that b = is equivalent to the 
u i,i+i = fJ-o- F° r other values of b /i is not const and thus the apparent mass is not const either. The distribution of 
the apparent mass on filled mass has maximum, with value which depends on b shown at Fig. {7\ 

One can see from this that there is no universal rescaled curve for different parameters. This fact was also mentioned 
in 0. Within the approach proposed this is the consequence of the existence of additional scale parameter X^. Such 
a scale becomes evident in the presence of overweight if /i is not constant but u is distributed within bulk. According 
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FIG. 7: Dependence the maximum of apparent mass on the fi changing gap in case of 80.83 g overweight presence. 

to the proposed distribution of the n, A^ depends on b and c, the gap in which /i changes and the steepness of fi 
saturation. 

The dependencies shown on Fig. must be parameterized using this parameter. In other words, this parameter 
splits Janssen's curve if silo is overloaded. 



We propose the description of the granular media at rest based on the introduction of the spatial distribution for 
contact coefficients n of dry friction. With the help of the simple static model we investigate the distribution of 
the weight in the silo. It is shown that in the case without overweight and homogeneous distribution of the friction 
coefficients the Janssen result is recovered. In a case of overweight we predict the maximum for the apparent mass as a 
function of a filled one, which is observed in experiments. It is important that the nature of such a maximum is related 
to the inhomogeneity in the spatial distribution of the dry friction coefficients. Such a distribution is formed due to 
the jamming of the grains in the upper layers which bear most of the overload. We put forward physical arguments 
which allow to obtain such distribution of the friction coefficients by taking into account the inhomogeneity of the 
pressure distribution with the height. Note the these result are obtained without any assumption about elasticity 
of the grains made in quasielastic approaches [|| . In addition within the proposed approach it is possible to explain 
the absence of simple rescaling law for the overshooting effect by the presence of two characteristic length. The first 
length is the Janssen length Aj and characterizes the distribution of pressure, the other one is the scale of spatial 
distribution for dry friction coefficient. 

In the continuous limit the development of the approach proposed implies the introduction of the tensor field for 
dry friction. Such a field becomes additional characteristic of a granular medium which is determined by both the 
geometrical and the stress configurations. The possibility of such a description is due to the fact that the interparticle 
interactions in the granular media are of short distance character It gives the grounds to expect that at least 
the static of granular medium should be described with the traditional framework of general elasticity theory with 
proper modifications. 



Authors thank to Dr. K. Shundyak for numerous discussions about the results. His critical remarks were also 
substantial for the final form of the paper. 
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